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GLOSSARY 


I.  INTRODUCTION 


A  "complete"  electromagnetic  pulse  (EMP)  hardness  statement  for  an 
aircraft  or  other  vulnerable  systems  requires  the  knowledge  of  both  the 
threat-induced  stress  or  its  bound  and  the  susceptibility  threshold  level  of 
every  pin  to  every  Line  Replaceable  Unit  (LRU).  The  induced  stresses  can  be 
obtained  by  measuring,  in  a  system-level  simulation  test,  the  currents  or  the 
voltages  on  the  wires  connected  to  LRUs  (and  then  using  the  appropriate 
extrapolation).  The  susceptibility  thresholds  can  be  obtained  by  performing 
onboard  and/or  laboratory  direct -drive  measurements.  However,  this  is  a 
prohibitive  effort  because  of  the  large  number  of  LRU  pins  contained  in  an 
ai rcraft. 

Attempts  have  been  made  continually  to  reduce  the  number  of  measurements 
required  for  a  reliable  hardness  assessment.  Since  wires  leading  to  LRU  pins 
are  usually  lumped  into  bundles,  one  such  attempt  was  to  find  relationships 
among  the  bounds  of  the  cable  bundle  current  and  the  bounds  of  the  individual 
wire  currents  by  using  a  statistical  approach  (Refs.  1,2).  However,  the 
samples  used  in  the  study  were  not  large  enough  and  the  conclusions  derived 
therefrom  were  not  convincing.  Furthermore,  the  bounding  approach  that  was 
used,  although  appropriate  for  pin  stress  estimates,  does  not  provide  pin 
susceptibility  thresholds. 

In  this  report,  a  deterministic,  as  well  as  a  more  general  approach,  is 
undertaken.  More  specifically,  analyses  are  performed  to  determine  how 
currents  or  voltages  on  individual  wires  of  a  bundle  can  be  established  with¬ 
out  measuring  each  wire  current  or  voltage.  It  is  expected  that  the  number  of 
required  measurements  will  strongly  depend  on  how  a  bundle  is  excited  and 
configured.  This  dependence  suggests  that  a  deterministic  approach  is  pre¬ 
ferred  over  a  statistical  one. 


In  Section  II,  a  detailed  analysis  is  given  as  to  how  the  individual  wire 
currents  can  be  determined  without  measuring  the  current  in  each  wire  .  The 
analysis  is  based  on  a  variational  technique.  Section  III  describes  how  the 
results  presented  in  Section  II  can  be  used  in  direct-drive  testing  in  general 
and  for  determining  pin  susceptibility  thresholds  in  particular.  This  type  of 
testing  can  be  utilized  for  developing  aircraft  electronic  component 
specifications  as  well  as  for  verifying  specification  requirements.  Sections 
II  and  III  conclude  with  simple  examples  that  illustrate  the  methodologies 
given.  Section  IV  is  devoted  to  a  detailed  study  of  the  currents  on  the  cable 
wires  and  how  the  excitation  mechanisms  and  cable  configurations  will  affect 
the  wire  responses.  Simple  and  useful  rules  for  selecting  the  optimal 
measurement  approach*  is  deduced  from  an  analysis  of  2-wire  bundles.  The 
variational  technique  introduced  in  Section  II  requires  information  on  the 
impedances  of  the  cable  bundle  at  the  pins.  Unfortunately,  the  existing  data 
base  does  not  provide  sufficient  information  to  verify  the  rules  deduced  from 
the  variational  technique.  Section  V  describes  new  experiments  for  the 
verification  and  improvement  of  the  rules.  The  experiments  are  designed  so 
that  tney  can  be  performed  either  in  a  laboratory  or  in  the  EMP  test  aircraft 
(EMPTAC).  The  conclusions  are  given  in  Section  VI. 


♦Optimality  here  is  defined  in  terms  of  the  number  of  measurements  required  to 
obtain  meaningful  estimates  of  each  wire  current. 


II.  DEDUCTION  OF  INDIVIDUAL  WIRE  CURRENTS  FROM  A  LIMITED  NUMBER 

OF  MEASUREMENTS 


When  an  aircraft  is  exposed  to  an  EMP,  high-intensity  transients  may  be 
induced  in  cable  bundles  inside  the  aircraft.  These  bundles  are  connected  to 
LRUs.  One  of  the  important  problems  relating  to  this  EMP  interaction  is  to 
find  the  intensities  of  the  induced  transients  at  the  pins  of  these  LRUs. 
These  transients  are  the  pin  currents  and  the  pin  voltages.  Obviously,  a 
procedure  for  quantifying  all  such  transients  based  solely  on  measurements  is 
not  practical.  Hence,  analytical  techniques  to  reduce  the  number  of  measure¬ 
ments  are  urgently  needed.  One  powerful  technique  is  that  of  a  variational 
principle,  which  requires  the  use  of  equivalent  circuits  at  the  LRUs.  The 
definition  of  the  equivalent  circuits  at  an  N-port  LRU  is  given  first,  and 
then  follows  the  application  of  a  variational  principle  to  the  problem  at 
hand. 


1.  DEFINITION  OF  AN  N-PORT  EQUIVALENT  CIRCUIT 

Let  a  cable  bundle  with  N  wires  be  connected  to  an  N-port  LRU.  This  cable 

guides  induced  transients  into  the  LRU  {Fig.  1).  The  driving  circuit  can  be 

replaced  by  an  equivalent  circuit,  namely  the  Thevenin  equivalent  circuit  with 

N  voltage  sources  which  can  be  denoted  by  a  vector  of  dimension  N  in  series 

with  an  N  x  N  source  impedance  matrix  (Zs  ).  The  load  impedance  (Z.  )  is 

n»m  Ln,m 

also  an  N  x  N  matrix  {Fig.  2).  Values  for  the  N  voltage  sources  contained  in 
this  circuit  can  be  obtained  by  disconnecting  all  the  wires  from  the  LRU  and 
then  measuring  all  the  voltages  at  the  disconnection  (Fig.  3a).  That  is 


<V  >  *  <V  Voc/ 

n  12 


•  ) 


oc. 


(1) 


where  (  stands  for  transpose  of  a  matrix,  and  Voc  is  the  open-circuit 


voltage  at  the  n-th  disconnection  while  all  the  wires  are  disconnected  from 


Figure  1.  An  N-port  LRU  with  an  N-wire  cable  bundle  connected  to  it. 


the  LRU.  The  first  step  taken  to  find  the  elements  of  the  source  impedance- 

matrix,  Z  ,  is  to  disconnect  all  the  wires  from  the  LRU,  apply  a  source  to 
sn,m 

the  m-th  wire,  and  then  measure  the  current  on  the  m-th  wire,  and  the  open- 

circuit  voltage  on  the  n-th  wire  while  the  rest  of  the  wires  remains  open 

(Fig.  3b).  With  this  procedure  Z,.  is  simply  given  by 

sn,m 


zs  '  V‘»  wBe"  *1  ■  h  ■ 

n  ,m 


•  •  •  •  "  1 


m-1  m+1 


=  ....  =  =  U  and  I  *  0. 

N  m 


(2) 


This  method  can  also  be  used  to  obtain  the  load  impedance  matrix  of  the  LRU 


Figure  3.  Experimental  setups  to  measure  the  open-circuit  voltages 
and  the  source  impedance  element  Z 

sn,m 


Another  equivalent  circuit,  namely  the  Norton  equivalent  circuit  with 
current  sources,  can  also  be  constructed.  The  current  source  vector  can  be 
found  by  shorting  out  all  the  wires  and  measuring  all  the  currents  flowing  on 
them,  regardless  of  the  nature  of  excitation  (Fig.  4a).  In  the  Norton  equiva¬ 
lent  circuit,  (Y  )  is  the  inverse  of  (Z  ),  i.e., 
sn,m  5n,m 

(Ys  )  *  (Z$  f1  (3) 

n,m  n  ,m 

Another  way  of  finding  the  element  Y  of  the  source  admittance  matrix  is  to 

sn,m 

disconnect  all  wires  from  the  LRU,  to  apply  a  source  to  the  m-th  wire,  and  to 


"'•  *  -  **- - 


10 


(a)  Short-circuit  currents  (b)  Source  admittance  elements  Y$ 

Figure  4.  Experimental  setups  to  measure  the  short-circuit  currents  n’m 
and  the  source  admittance  elements  Y 

sn,m 

measure  the  voltage  on  the  m-th  wire  and  the  short-circuit  current  on  the  n-th 

wire  while  the  rest  of  the  wires  remains  shorted  (Fig.  4b).  Then 

Y  is  simply  given  by 
sn,m 


Y 

sn  ,m 


I  /V  when  V,  =  V, 
n  m  IE 


V  =  0  and  V  *  0  (4) 
N  m  ' 


From  the  Thevenin  equivalent  circuits,  write 


(I,  )  *  (U,  )  +  Us  )fl 

n  n,m  n,m 


) 


-1 


(5) 


(V,  )  -  (Z.  )  •  ((Z.  )  +  (Ze  )) 


(V  ) 


(6) 


Similarly,  from  the  Norton  equivalent  circuit,  write 


'V  ■  •  “sc  >  <7> 

n  n,m  n  ,m  m 

(IL  )  ■  (VL  )  •  ((Yl  )  +  Y  ))_1  .  (I  )  (8) 

n  n,m  m,p  m,p  p 

Since  these  two  equivalent  circuits  should  give  the  same  load  responses,  it  is 

obtained,  for  (Z.  )  =  (0  ) 

Ln,i» 

(I,  )  ■  (I„  )  ■  (2,  I'1  •  (*  )  ■  (V,  )  •  (»„  )  (9) 


The  load  responses,  (V,  )  and  (I.  ),  are  the  quantities  of  interest. 

Ln  Ln 

However,  it  is  very  difficult,  if  not  practically  impossible,  to  find  all  the 

elements  of  (IL  )  and  (V^  )  for  the  LRU.  This  is  due  to  the  fact  that  to 
n  n 

find  (I,  )  or  (V,  )  for  all  excitations  necessary  for  a  reliable  hardness 
n  n 

assessment  requires  a  large  number  of  measurements.  It  is  therefore  desirable 
to  quantify  pin  currents  and  pin  voltages  with  fewer  number  of  measurements. 

To  make  this  possible,  one  needs  to  find  the  relationships  among  the  fewer 
measured  quantities  and  all  the  individual  pin  responses.  The  fewer  measure¬ 
ments  should  be  the  ones  that  are  easy  to  measure  and  expected  to  contain  more 
information.  Some  of  the  candidate  quantities  for  measurements  are  the  cable 
bundle  current  IB  and  the  currents  on  wires  that  are  expected  to  be  strongly 
excited.  In  the  following,  the  relationships  will  be  developed  when  only  the 
cable  bundle  current  is  measured,  or  when  only  one  wire  current  is  measured, 
or  when  one  wire  current  and  the  cable  bundle  current  are  measured.  The 
relationships  can  then  easily  be  extended  to  cases  with  more  number  of  mea¬ 


surements. 


2.  ONLY  Ig  MEASURED 


It  is  known  that  currents  and  voltages  in  a  circuit  distribute  themselves 
in  such  a  way  that  the  energy  is  stationary  with  respect  to  currents 
(Ref.  3).  In  addition  to  satisfying  the  stationary  condition,  constraints 
must  be  imposed  for  determining  all  the  currents/ voltages.  The  only 
constraint  to  be  considered  in  this  subsection  is  that  the  sum  of  the  individ¬ 
ual  wire  currents  be  equal  to  the  measured  bulk  cable  current  Ig.  The  method 
of  Lagrange  multipliers  is  now  used  to  find  the  relationship  between  Ig  and 
the  individual  wire  currents  (Refs.  4,  5,  6).  The  energy  function  of  the 
Thevenin  equivalent  circuit  is  given  by 


p  *  (IL  )  •  (Zt  )  .  (IL  ) 
n  n,m  m 

(10) 

where  (Zt  )  *  (ZL  )  +  (Z$  )  which  are  generally  complex  matrices. 

n,m  n,m  n,m 

Rewrite  P  in  terms  of  the  elements  of  the  matrices.  That  is 

N  N  * 

p  •  l  l  li 

n«l  m-1  n  n,m  m 

(ID 

The  constraint  is 

N 

=  ill  Ilm 

(12) 

with  I.  being  the  m-th  pin  current, 
m 

Substituting 

I,  *  a  +  jb 
m 

(13) 

*  LBr  +  jIBim 

into  Equations  11  and  12  obtains 

(14) 

N  N 

l  l 

n*l  m*l 


(a 


JVZt 


n  .m 


(am+ 


JV 


(15) 


'*■ '  Jl a" 

N 

*8im  bm 

m=l 


Now,  we  must  find  the  stationary  points  of  the  expression 

'  *  j-i  1 <a"  ■  »(a" + jt>")  * x^'  i,  a">! 

n=i  m=l  n,m  m=l 

"  ,bm>  <18> 

m=l 

where  A ^  and  Ap  are  Lagrange  multipliers.  Accordingly,  we  should  have 

U-  ■  0,  n  -  1,  2,  ....  N  (19) 

3an 

|f-  *  0,  n  »  1,  2,  ....  N  (20) 

Consequently,  we  get 

Jl  Zt„  <*„  *  *  j.  <*.  *  *  ‘l  *  0  (21> 

m=l  n,m  m=l  m,n 

j,  \  mlam  *  Jbm>  *  j.  <»„  *  Jb„>Zt„,  ’  x2  ’  0  <22> 

m=i  n,m  m= l  m,n 

Assuming  that  (Z  )  is  a  symmetric  matrix,  one  can  rewrite  Equations  21  and 


22  as 


2<Zt  >  *  <am>  -  *1<V 

n  ,m 

2<Zt  >  *  <V  -  W 

n  ,m 


where 


Um)T  *  ^al»  a2*  *  *  *  *3n)  *  (bm)T  *  (bi»  b2»  *  *  *  »bfl) »  and  Un)1  *  U»  l>  •••»!) 


(a  )  and  (b  )  can  thus  be  written  as  follows 
m  m 


f  *-'**-■  iw\  In  i/i  fc*  hi fci 


V.VA  V  V  V  N.  V*  “A  V  ^  ^  "7,  -  .  •  «  H  -  M"-  <V 


,vr:.*-- .Vi.*-. r— •x.»  ._*  h.-.  - 1 


(aj4h(\  >  •  <V 

m,n 


^>  -7‘2<Vt  »  •  <‘nJ  <Z6> 

m,n 

where  (Y  )  *  (Z^  )_1.  Substituting  Equations  25  and  26  into  Equation  16 

tn,m  ^n.m 


and  17,  and  solving  for  and  we  find 


A  <V •  <vt  >-<V 

n,m 

2IU 

X  n  =  - r - — -  (28) 

(!„)  •  (’t  )  ■  <V 

n  ,m 

Consequently,  by  substituting  into  Equations  25  and  26  (am)  and  (bm)  are  found 


to  be 


(V  *  *Br 


(Yt  )  •  <y 

n,m 


<V  *  <V  >  *  (V 

n  ,m 


<bm>  *  ‘aim 


<V  J  '  l1"’ 

n  ,in 

(‘n)'  '  <V  >  '  (V 

n  ,m 


Finally,  we  have  arrived  at  the  relationship  between  lg  and  (IL  ), 

<V 

(I,  )  ■  I» - 1— 1 -  (31 

”  <V  •  <  Yt  >  '  (V 

n  ,m 

where 

(Yt  )  *  (Z  rl  •  <(Z  )  ♦  (Z  ))_1 

m,n  m,n  m,n  m,n 

IR  is  the  measured  bulk  cable  current  and  Ii  is  the  current  on  wire  m.  As 
0 

can  be  seen,  from  a  knowledge  of  Ig,  (Z$  )  and  (Z^  ),  (IL  )  can  be 

n,m  n,m  m 

obtained.  Using  the  Thevenin  equivalent  circuit,  we  can  write 


(32) 


(VQr  )  -  (Z  )  •  (I,  )  •  IB  - , - - - 

°Cn  Vm  Lm  B  (ln)T  .  (Yt  )  •  (1J 


n,m 


Therefore,  (IL  )  obtained  in  Equation  31,  will  be  the  exact  current  vector,  if 
m 

the  open-circuit  voltages  on  all  the  wires  are  the  same.  Equations  31  and  32 


hold  true  in  the  time  domain  when  (Z,.  )  and  (ZL  )  are  resistive  for  all 


'm.n 


'm,n 


frequencies. 


3.  ONLY  ONE  INDIVIDUAL  WIKE  CURRENT  MEASURED 

Another  type  of  constraint  that  can  be  imposed  is  that  the  current  of  the 
most  strongly  excited  individual  wire,  for  example,  Ik  be  given  from  measure¬ 
ment  .  That  i s 

(33) 

From  Equations  13  and  33,  we  obtain 

(34) 

(35) 

where  Ikr  and  Ikim  are  real  and  imaginary  parts  of  Ik,  respectively.  Now,  we 
must  find  the  stationary  points  of  the  expression 


:L  *  Xk 
Lk 


ak  *  Xkr 


bk  *  *kim 


N  N 


l  {an  -  jbn>Zt  ffl(*m  +  jbm>  +  h^kr  *  ak>  +  ^kim  ‘  V  (36> 
n=l  m=l  n,m 


where  Xj  and  X2  are  Lagrange  Multipliers.  Accordingly  we  should  have,  from 
Equations  19  ana  20, 


1  Zt_  *  JbJ  *  I  l*.  -  Jb„>Zt_  .  *  MnK  *  0  <37> 

m=l 


n=l  n,m 


m,n 


•J  J,  Zt  <am  *  Jbm>  *Jj,  (am  *  jbn,>Zt  .  *  Vnk  ‘  0  <38> 


m*l  n,m 


m=l 


m,n 
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Assuming  that  (Z  )  is  a  symmetric  matrix,  one  can  write  Equations  37  and  38 
m,n 
as 


2<Zt  >  •  <*.>  *  xl‘nk 
n  ,m 

(39) 

2<zt  >  •  (6m>  ■  x2s„k 

n  ,m 

(40) 

a  and  b  can  thus  be  written  as  follows 
m  m 

3  *  i  X,Y 

m  <•  1  t  v 

m,K 

(41) 

bm  =  7X2Yt  k 
m,k 

(42) 

Substituting  Equations  41  and  42  into  Equations  34  and  35,  and  solving 


for  Xj  and  x 2  we  find 


21 


kr 


Lk,k 


21 


X2  = 


kim 


Therefore,  a  and  b_  are 
m  m 


k,k 

b 

n 

Y. 


t  -  t  nitk 
am  "  *kr  T 


k,k 


bm  =  ^im  Y 


'm,k 


uk  ,k 


Finally,  we  arrive  at  the  following  important  relationship: 

I. 


I, 


L  *  lk  Y 
m  t 


k,k 


Using  the  Thevenin  equivalent  circuit,  we  get 


(43) 

(44) 


(45) 


(4b) 


(47) 


(yoc  >  *  <Zt  )  •  <*L  ) 
n  n,m  m 


Hence,  (1^  )  obtained  from  the  above  method  will  be  accurate  if  the  open- 
m 

circuit  voltages  of  all  wires  except  the  k-th  wire  are  zero.  Equations  47  and 

48  hold  true  in  the  time  domain  when  ( Z.  )  and  (Z«.  )  are  resistive  for  all 

Tn,n  sm,n 

frequencies . 

4.  ONLY  BULK-CABLE  CURRENT  ANO  CURRENT  ON  ONE  STRONGLY  EXCITED  WIRE  MEASURED 
Using  an  approach  similar  to  the  ones  used  in  the  previous  cases,  one 


obtains. 


I,  *  h  “ 

Ln  B  Y 


N  N 

Yt  I  Yt  •  Yt  l  Yt 

Lk,k  m=l  n,m  n,k  m=l  ’'k.m 

__  _  j_  |f 


II  My  H  y 

tk-k  ml  p=l  '"-P  ml  tm-k 


V 

l  \  m 
m-1  k’m 


N  N  N  N 

Yt  I  l  Y  -  l  yt  l  Yt 

T  n,k  m=l  p=l  xm,p  m*l  m,k  m=l  n,m 

k  “  FI  IT~  R  "  R  Z  1  1 

tk-k  Jl  p-l  'm-P  ml  tn,-k  ml  tk-" 


n  y 

l.  \ 


where  Ig  and  Ik  are  the  measured  bulk-cable  and  the  k-th  wire  current,  respec¬ 
tively  (see  Appendix  A). 

The  above  method  can  be  applied  to  any  other  combination  of  current  mea¬ 
surements  as  constraints  for  obtaining  relationship  among  wire  currents. 

5.  EXAMPLES. 

Figure  5  depicts  a  two-wire  configuration.  The  wires  all  have  the  same 

length  ,  but  they  do  not  necessarily  have  the  same  displacement  above  the 

ground  plane.  The  characteristic  impedance  matrix  (Z„  )  provides 

'"n.m 
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information  aoout  the  line  configuration.  These  wires  are  terminated  at  botn 

ends  with  the  load  and  terminating  impedance  matrices,  (Z>  )  and  (Zn  ), 

'-n.m  °n,m 

respectively.  These  wires  are  assumed  to  be  excited  by  voltage  and  current 

sources  denoted  by  (V  )  and  ( I  ).  Computer  codes  have  been  developed  to 

n  n 

calculate  currents  and  voltages  at  any  point  on  these  wires  for  arbitrary 
sources,  load,  terminating,  and  characteristic  impedances. 

Let  the  sources  take  the  following  form: 


Vo  =  Vo  exp(‘  at)sin<2 
n  n 


IQ  (t)  =  IQ  exp(-  ot)sin(2  irfQt)  (bl) 

n  n 

where  f  =  Z  MHz  ,  a  =  5  x  106  s"1,  V  =  1  V  and  I_  =  1  mA.  The  chosen 
0  °n  un 

values  for  f  and  a  are  typical  for  the  fundamental  resonant  mode  of  most 
aircraft  (Ref.  7).  Using  the  computer  code,  Ij  and  I2  are  found  at 
xm  *  s.  .  These  exact  quantities  are  now  compared  with  those  obtained  from  the 
variational  technique  described  previously.  First,  take  the  case  where  only 
Ig  is  measured,  whose  value  is  just  1^  +  I2.  Using  this  Ig  value  in  Equation 
31  obtains  estimates  of  the  wire  currents,  1^  and  1^.  Next,  consider  the  case 
where  only  the  current  on  wire  #1  is  measured.  Using  this  current  value  as  Ik 
in  Equation  47,  obtains  estimates  of  wire  currents, 

I"  and  I". 

Figures  6  and  7  present  the  plots  of  these  currents  for  two  different 

types  of  sources.  Since  the  terminating  impedance,  (Z  )  at  x  s  0  is  equal 

°n,m 

to  the  characteristic  impedance  (Z  )  ,  the  open-circuit  voltages  of  the  two 

cn,m 

wires  are  about  the  same.  Consequently,  1^  and  l£  obtained  from  Ig,  alone  are 
very  much  similar  to  and  I2  as  they  should.  However,  is  at  variance 
with  I2.  To  quantify  the  accuracy  associated  with  the  currents  estimated  from 
only  a  limited  number  of  measurements,  a  curve  is  introduced,  as  shown  in 


m 


exp(- at)sin(2ir fQt)  mA 
*  5  x  106  s’1 


a 


Figure  8,  in  which  the  horizontal  axis  indicates  che  number  of  measurements  in 
a  bundle  of  N  wires  and  the  vertical  axis  provides  a  measure  of  accuracy  on 
the  individual  wire  currents.  The  vertical  axis  extends  from  0  to  100 


percent.  A  vertical  value  of  100  percent  indicates  that  there  is  no 
difference  in  the  estimated  value  and  the  actual  value.  One  accuracy  gauge  to 


be  used  for  the  vertical  axis  is 


Accuracy  =  a„  wi™  £■  cot' 


I  max 1 1 


actual  1 
max  1 1 


max  1 1 


calc. 1 


actual 


x  100%  (52) 


Figure  9  depicts  the  accuracy  in  terms  of  the  relative  error  which  is  the 
argument  of  tne  inverse  cotangent  function.  As  can  be  seen,  the  relative 
error  of  100%,  i.e.,  |Ica1c  I  *  2| ^actual  ^  corresP°nds  t0  an  accuracy  of  50 
percent. 

It  is  expected  that  the  accuracy  curve  defined  by  Equation  52  be  monoto- 
nically  increasing  with  the  number  of  measurements.  In  the  two-wire 
transmission-line  example  the  accuracy  is  100  percent  if  the  number  of 
measurement  is  2.  However,  if  the  number  of  measurement  is  1,  e.g.,  only  an 
individual  wire  current  is  measured,  then  a  relative  error  of  62  percent  is 
obtained  from  Figure  6b,  which  corresponds  to  accuracy  =  64  percent. 

Figures  10  and  11  show  the  cases  where  only  one  wire  is  excited  and  the 


terminating  impedance  (Z  )  is  equal  to  the  characteristic  impedance 

°n,m 

(Z  ).  In  these  cases,  the  open-circuit  voltage  of  the  excited  wire  is 
cn,m 

nonzero  and  dominant  over  the  others.  Therefore,  I and  are  very  much 
similar  to  Ij  and  I 2»  as  can  be  seen  from  Figures  10  and  11.  However, 

1^  and  I£  are  not  equal  to  1^  and  l£.  The  corresponding  relative  error  and 
accuracy  for  Ij^  in  Figure  11  are  75  and  59  percent,  respectively. 

Figures  12  and  13  show  the  cases  where  the  terminating  impedance 


(Z  )  at  x  =  0  is  not  equal  to  the  characteristic  impedance  (Zc  ).  As  can 
^  jin  n  jdi 

be  seen  from  Figure  12  where  the  two  wires  are  excited,  Ij  and  I), 
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NUMBER  OF  MEASUREMENTS 


Figure  8.  The  accuracy  measure  versus  the  number  of  measurements 
for  two-wire  transmission  line  model. 

obtained  from  IB  have  a  smaller  maximum  relative  error.  On  the  other  hand, 
where  one  wire  is  excited  (Fig.  13),  I”  and  I"  have  a  smaller  maximum  relative 
error. 

Additional  curves  presenting  the  comparison  of  exact  wire  currents  with 
calculated  ones  are  given  in  Appendix  B. 

In  an  N-wire  cable,  depending  on  the  measurement  sequence,  the  rate  of 
increase  of  the  accuracy  curve  will  be  different.  It  is  possible  to  obtain  an 
accuracy  greater  than  a  certain  value  from  a  number  of  measurements  much  less 
that  N. 

The  following  conclusions  can  be  made  from  the  above  observations 
concerning  multi wire  cables: 

•  Wire  currents  can  be  calculated  from  a  limited  number  of  measurements. 


k 


RELATIVE  ERROR  (%) 


Figure  9.  The  accuracy  measure  versus  the  relative  error. 


•  If  the  wires  are  evenly  excited,  wire  currents  calculated  from  a 
knowledge  of  the  bulk-cable  current  have  a  higher  accuracy  than  that 
obtained  from  other  kinds  of  current  measurements. 

•  If  the  excitation  of  one  individual  wire  is  dominant  over  the  others, 
wire  currents  calculated  from  a  knowledge  of  the  current  of  that 
specific  wire  have  a  higher  accuracy  than  that  obtained  from  other 
kinds  of  current  measurements. 

The  variational  technique  introduced  in  this  section  can  be  applied  to 
obtain  useful  information  on  the  wire  currents  from  a  limited  number  of  mea¬ 
surements.  Although,  the  technique  requires  information  on  the  source  as  well 
as  the  load  impedance  of  the  cable  bundle  at  the  point  of  interest,  the  imped¬ 
ance  measurements  can  be  performed  without  the  need  of  bringing  the  aircraft 
to  the  EMP  simulator  site.  Another  point  is  that  impedance  measurement  is 


•  *  *  *  •  */  *  *  v  >v  v  .*  »•  -* 
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II 

V- 


(0'K,M«: 


9  48.3 

3  317.9, 


(!o  )  =  (o)  6XP^"  at)sin(2»f0t) 


a  =  5  x  106  s"1 

f  *  2  MHz 
0 


Jjexp(-  at)sin(2» fQt)  ^  Q 

a  ■  5  x  106  s'1  /  \  _\ 

(zv.)'(°  50)°  • 


I  i  nc.  k  jj  s  i 
)  Time  history  of  1^.  I^.  and 

act  wire  currents  (Ij,  1^)  and  estimated  wire  curr 
1»  ij;  I'j,  I^)  using  two  variational  principle  tech 
en  L  =  75  m,  x$  =  18.75  m,  xffl  =  75  m,  Z^j  =  100  ft 
12  =  00 ,  and  Z„-2  =  400  ft. 


III.  DIRECT-DRIVE  TESTING 


As  stated  in  the  Introduction  a  “complete"  EMP  hardness  statement  for  an 
aircraft  requires  a  knowledge  of  threat-induced  stress  or  its  bound  and 
susceptibility  threshold  for  each  pin  of  each  LRU.  Direct-drive  testing  is 
typically  the  method  by  which  susceptibility  thresholds  can  be  obtained.  This 
type  of  test  can  also  be  used  to  verify  if  components  and  subsystems  meet 
their  required  specifications.  The  results  obtained  in  Section  II  are  applied 
in  this  section  to  the  development  of  two  test  procedures  for  direct-drive 
testing.  One  procedure  is  the  excitation  of  the  total  bundle  by  a  direct- 
drive.  The  other  is  direct  driving  of  a  single  wire  in  the  bundle. 


1.  DIRECT  DRIVING  OF  CABLE  BUNDLE 

Consider  a  cable  bundle  with  N  which  are  terminated  at  one  end  at  an 
LRU.  Figure  14  shows  an  inverse  current  probe  exciting  the  entire  bundle.  It 
is  observed  that  each  wire  contained  in  the  bundle  is  excited  by  the  same 
amount  of  source  magnetic  flux.  Consequently,  the  open-circuit  voltages,  Voc, 
induced  by  the  direct -drive  excitation  are  the  same  (Ref.  9).  Section  II 
shows  that  when  the  VQCs  of  the  equivalent  circuits  were  about  the  same,  the 
bulk  cable  current  Ig  would  be  very  informative  in  evaluating  individual  wire 
currents.  In  fact,  if  Ig  is  measured,  then  the  individual  wire  currents  can 
be  found  from  Equation  31,  viz.. 


m,n 


*1  =  - T— 

■»  (1J  .  (  )  •  (1 J 


(53) 


n,m 


m' 


Equation  53  can  be  used  in  verification  of  LRU  specifications.  Consider 
an  LRU  in  a  laboratory.  This  LRU  can  be  connected  by  an  N-wire  cable  bundle 
to  another  box  with  an  impedance  matrix  (Zn  ).  It  is  assumed  that  the 


'n,m 


■1 /•'V.-V-O  .. 


"V  • 
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Figure  14.  Total  cable  bundle  excitation  using  an  inverse  current  probe. 


impedance  matrix  (Z  )  can  be  changed.  To  verify  the  pin  specifications  of 
°n,m 

a  given  LRU,  1D  and  the  elements  of  the  impedance  matrix  Z  can  be 

o  . 

n  ,m 

adjusted  until  the  current  or  voltage  of  each  pin  becomes  identical  to  the  pin 

specification.  Therefore,  by  using  an  inverse  current  probe  and  a  proper 

impedance  matrix  (Zn  ),  it  is  possible  to  verify  a  pin  specification. 

°n,m 

It  is  interesting  to  note  that  Equation  53  can  also  be  used  in  defining 
the  specification  of  an  LRU  in  terms  of  Ig.  To  do  so,  using  Equation  53,  one 
can  obtain  the  specification  of  LRU  from  that  of  the  sensitive  pin.  There¬ 
fore,  specification  can  be  verified  by  direct -driving  the  cable  bundle  using 
an  inverse  current  probe. 


2.  DIRECT  DRIVING  OF  A  SINGLE  WIRE 


In  this  example,  only  one  wire  of  a  cable  bundle  is  directly  driven. 
Therefore,  the  VQC  of  this  wire  is  usually  dominant  over  the  open-circuit 
voltages  of  the  other  wires.  Based  on  the  results  obtained  in  Section  II,  the 
current  in  the  excited  wire  can  provide  information  about  the  other  wire 
currents. 

If  the  current  in  the  excited  wire,  which  is  denoted  by  1^,  is  measured 
the  other  wire  currents  can  be  computed  from 


I,  =  I 


k  T 


m,k 


(54) 


*k,k 


This  approach  can  also  reduce  the  number  of  measurements  required  to  quantify 
pin  susceptibility  thresholds.  To  further  illustrate  this  point,  consider  an 
example  of  a  four-port  LRU  connected  to  a  driving  circuit  through  a  four-wire 
cable  bundle.  The  total  admittance  of  the  LRU  is  assumed  to  be 


<Yt 

Lm,n 


) 


1.82 

-1.55 

1.75 

-0.12  \ 

X 

*— • 

O 

1 

CO 

§ 

0 

-1.55 

2.31 

-1.92 

2.5 

1.75 

-1.92 

1.65 

-1.3  j 

-0.12 

2.5 

-1.3 

2.01  / 

(55) 


As  can  be  easily  seen  only  the  third  column  or  row  has  two  off-diagonal  ele¬ 
ments  (1.75,-1.92)  whose  absolute  values  are  greater  than  the  diagonal  term 
(1.65).  This  implies  that  if  wire  No.  3  is  excited,  according  to  Equation  54, 
wires  Nos.  1  and  2  will  reach  their  maximum  allowable  values  before  wire  No.  3 
will.  If  the  susceptibility  thresholds  of  all  the  three  wires  are,  say,  100 
mA,  then  one  can  calculate  from  Equations  54  and  55  that  for  I|_^  *  100  mA,  I3 
should  be  94.2  mA,  which  is  less  than  100  mA.  Similarly  for  I  =  100  mA,  I3 
turns  out  to  be  85.9  mA,  which  is  also  less  than  100  mA.  Therefore,  without 
changing  the  experimental  setups,  i.e.,  direct  driving  wire  No.  3  alone,  the 
susceptibility  thresholds  of  the  first  two  pins  can  be  established. 


IV.  EFFECT  OF  VARIOUS  PARAMETERS  ON  COMMON-  AND  DIFFERENTIAL-MODE  CURRENTS 


To  select  an  optimal  measurement  approach,  the  effects  of  different  param¬ 
eters  on  wire  currents  must  be  studied.  In  this  section,  the  effects  of 
various  parameters  (such  as  the  excitation  mechanism,  bundle  configuration, 
impedances,  source  and  measurement  locations,  etc.)  on  wire  currents  are 
studied.  Before  presenting  the  results,  the  concept  of  common  and  differen¬ 
tial  modes  is  introduced. 


1.  DEFINITION  OF  COMMON  AND  DIFFERENTIAL  MODES 

As  mentioned  earlier,  the  sum  of  wire  currents  is  referred  to  as  the  bulk- 
cable  current  Ig  and  it  is  relatively  easy  to  measure.  However,  Ig  gives  only 
partial  information  about  the  individual  wire  currents.  Additional 
information  is  needed  to  obtain  a  complete  description  of  the  wire  currents. 
Tnis  additional  information  can  be  other  combinations  of  wire  currents.  For 
example,  for  a  two-wire  cable,  Ig  is  the  sum  of  Ij  and  Ig;  i.e.,  Ig  =  Ii  +  Ig, 
and  the  other  combination  is  the  difference  between  1^  and  I^,  i.e.. 

Ip  *  I ^  -  I^.  Ig  and  Ip  may  be  referred  to  as  the  common  and  differential 
modes.  For  the  three-wire  cable,  the  common  and  differential  modes  are  writ¬ 
ten  as: 


*  *1  +  *2  +  *3 
*D1  =  *1  +  l2  "  *3 
l02  =  l\  ‘  lZ  +  *3 


Transmission  line  models  are  used  in  the  analysis  that  follows  to  study 
the  effects  of  excitations,  configurations,  etc.,  on  bulk  and  di fferenti al - 
mode  currents. 
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2.  EFFECT  OF  LOAD  IMPEDANCE 

Consider  a  two-wire  line  with  the  terminating  impedance  matched  to  the 
characteristic  impedance.  Ig  and  Iq  are  calculated  for  each  value  of  the  load 
impedance.  Figures  15  through  17  present  these  results.  The  following 
conclusions  can  be  drawn  from  these  figures: 

•  For  (two)  like  pins  (i.e.  the  pin-to-ground  impedances  of  the  two  pins 
are  the  same),  Iq  is  sensitive  to  pin-pin  impedance,  whereas  Ig  is 
not. 

•  For  (two)  unlike  pins  (i.e.,  the  pin-to-ground  impedances  of  pins  are 
not  the  same),  Ig  and  Iq  do  not,  in  general,  depend  appreciably  on  the 
pin-pin  impedance. 

More  examples  of  Ig  and  Iq  for  different  load  impedances  are  presented  in 
Appendix  B. 

3.  EFFECT  OF  SOURCE  TYPE 

In  this  case,  Ig  and  Iq  are  calculated  for  different  sources:  voltage  or 
current.  Figure  18  presents  Ig  and  Iq  for  different  sources.  The  following 
conclusions  can  be  drawn  from  Figure  18: 

•  Regardless  of  the  load  and  source  impedances,  the  Ib/Iq  interaction  is 
independent  of  the  type  of  source,  i.e.,  a  voltage  or  current 
source.  This  leads  to  the  following  corollaries: 

-  If  Ig  (or  Iq)  is  large  compared  to  Iq  (or  Ig)  with  one  type  of 
source,  it  is  also  large  for  the  other  type  of  source. 

-  If  one  type  of  source  gives  rise  to  zero  Ig  (or  Iq),  so  does  the 

other. 

-  For  direct -drive  tests,  either  the  voltage  or  current  source  will 
produce  the  same  wire  current  distribution. 

Additional  figures  of  Ig  and  Iq  for  sources  are  given  in  Appendix  B. 


EFFECT  OF  LORD  I MPEDRNCE 
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4.  EFFECT  OF  SOURCE  POSITION 

For  this  case,  sources  are  located  at  different  positions  and  Ig  and  Ig 
are  calculated  for  each  position.  Figure  19  shows  Ig  and  Ig  for  different 
source  positions.  The  following  conclusions  can  Oe  drawn  from  this  figure. 

•  Regardless  of  load  impedance,  changing  the  position  of  the  source 
merely  results  in  a  shift  of  the  waveforms  of  Ig  and  Ig.  Therefore, 
the  Ig/ Ig  interaction  is,  indepedent  of  source  positions.  For  direct 
drive  tests,  the  location  of  the  drive  can  be  safely  assumed  to  have 
no  effect  on  the  wire  current  distribution  at  the  load. 

There  are  additional  figures  in  Appendix  B  showing  Ig  and  lg  for  different 
source  positions. 

5.  EFFECT  OF  NUMBER  OF  EXCITED  WIRES 

Figures  20  through  22  present  the  results  for  Ig  and  Ig  for  different 
numbers  of  locally  excited  wires.  The  conclusions  drawn  are  as  follows: 

•  For  (two)  like  pins,  changing  the  number  of  excited  wires  from  one  to 
two  generally  causes  Ig  to  increase  and  Ig  to  drop  to  zero  regardless 
of  the  load  and  source  impedances. 

•  For  (two)  unlike  pins,  changing  the  number  of  excited  wires  from  one 
to  two  causes  Ig  to  increase  and  Ig  to  decrease  regardless  of  load  and 
source  impedances. 

There  are  additional  figures  in  Appendix  B  presenting  Ig  and  Ig  for  different 
numbers  of  excited  wires. 

6.  EFFECT  OF  LINE  CONFIGURATION 

The  effect  of  line  configuration  parameters  on  Ig  and  Ig  can  be  studied  by 
varying  the  elements  in  the  characteristic  impedance  matrix.  Figure  23 
presents  Ig  and  Ig  for  different  characteristic  impedances.  The  following 
important  conclusion  is  reached  from  study  of  these  figures: 

•  The  line  parameters  do  not  have  any  significant  effect  on  Ig  and  Ig. 
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An  additional  figure  of  1B  and  ID  for  different  characteristic  impedances  is 
given  in  Appendix  B. 

7.  EFFECT  OF  MEASUREMENT  LOCATION 

In  this  case,  Ig  and  Ig  are  calculated  at  different  locations  of  the 
line.  Figure  24  shows  these  results.  The  following  conclusion  is  drawn  from 
Figure  24: 

•  Regardless  of  the  load  impedance,  changing  the  measurement  point 
produces  only  a  shift  of  the  waveforms  of  IB  and  I,,. 


V.  SUGGESTED  EXPERIMENTS 


The  variational  technique  introduced  in  Section  II  requires  the  informa¬ 
tion  on  the  impedances  of  the  cable  bundle  at  the  pins.  Since  the  existing 
data  base  does  not  provide  such  information,  new  experiments  are  necessary  for 
the  verification  of  the  rules  established  in  previous  sections.  By  using  an 
accuracy  gauge,  such  as  that  defined  in  Equation  52,  the  verification  of  the 
rules  can  be  accompl i shed.  For  example,  take  an  accuracy  measurement  of 
greater  than  80%  to  be  a  satisfactory  criterion.  Then,  if  the  results 
obtained  from  the  rules  give  such  an  accuracy  measurement,  the  rules  can  be 
considered  to  be  valid.  Two  examplary  multiwire  experiments  are  suggested  to 
verify  the  following  rules: 

•  If  the  wires  are  evenly  excited,  wire  currents  calculated  from  a 
knowledge  of  the  bulk -cable  current  alone  can  have  a  satisfactory 
accuracy. 

•  If  the  excitation  of  one  individual  wire  is  dominant  over  the  others, 
wire  currents  calculated  from  a  knowledge  of  the  current  of  that 
specific  wire  can  have  a  satisfactory  accuracy. 

If  the  rules  do  not  provide  satisfactory  results  when  compared  with  the  exper¬ 
imental  results,  it  generally  means  that  additional  constraints  are  needed. 

To  demonstrate  this  point,  a  third  experiment  is  suggested  which  allows  for 
two  contraints. 

1.  EXPERIMENT  TU  VERIFY  THE  RULE  OF  USING  THE  VARIATIONAL  TECHNIQUE  AND  ONLY 

IB  MEASUREMENT 

In  an  aircraft  under  a  simulator  environment,  select  a  cable  bundle  whose 
wires  all  lead  to  points  of  entry  (POE)  of  similar  penetration  strength.  Then 
measure  IB  and  load  and  source  impedances  and  use  the  formula  obtained  in 
Section  II  (Equation  31)  to  calculate  the  wire  currents.  The  calculated  wire 
currents  are  then  compared  with  the  measured  ones.  The  degree  of  accuracy  can 
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then  be  calculated  and  the  verification  of  the  rule  for  the  selected  cable  can 
then  be  accomplished.  If  it  is  impossible  to  locate  such  a  cable  inside  an 
aircraft,  a  laboratory  setup  should  be  used  as  an  alternative. 

2.  EXPERIMENT  TO  VERIFY  THE  RULE  OF  USING  THE  VARIATIONAL  TECHNIQUE  AND  ONLY 

A  MEASUREMENT  OF  ONE  INDIVIDUAL  WIRE  CURRENT 

In  an  aircraft  under  a  simulator  environment,  select  a  cable  bundle  in 
which  only  one  wire  is  strongly  excited  by  a  POE.  Then  measure  the  current  on 
that  wire  and  load  and  source  impedances  and  use  Equation  47  to  calculate  each 
individual  wire  current.  The  calculated  wire  currents  are  then  compared  with 
the  measured  ones  to  obtain  the  accuracy  measure  (Equation  52).  The 
verification  of  the  rule  can  then  be  accomplished.  Again,  if  such  a  cable 
cannot  be  found  in  an  aircraft,  a  laboratory  setup  must  be  used  as  an 
alternative. 

3.  EXPERIMENT  TO  VERIFY  THE  RULE  OF  USING  THE  VARIATIONAL  TECHNIQUE  AND  TWO 

WIRE  CURRENT  MEASUREMENTS 

For  the  cable  bundles  selected  in  the  above  two  experiments,  it  is  possi¬ 
ble  that  the  accuracy  measures  obtained  are  too  low  to  be  satisfactory.  Under 
this  situation,  additional  measurements  should  be  made  (in  experiment  1, 
current  measurement  of  a  relatively  strongly  excited  wire;  in  experiment  2, 
measurement  of  I8)  and  Equation  49  applied  to  calculate  the  individual  wire 
currents.  These  calculated  currents  are  then  compared  with  the  measured 
currents  for  an  accuracy  measure. 

4.  EXTENDED  EXPERIMENT  TO  IDENTIFY  BEST  MEASUREMENT  COMBINATION 

The  previous  experiments  can  be  extended  to  determine  the  best  combination 


of  current  measurements  which  provides  a  satisfactory  accuracy  measure.  To  do 
so,  select  a  cable  bundle  in  the  aircraft.  For  the  selected  cable  bundle, 
first  determine  whether  the  rules  identified  in  the  previous  experiments  can 


De  used  and  perform  the  correspondi ny  experiments.  If  the  accuracy  measure 
obtained  therefrom  is  satisfactory,  the  search  for  the  best  measurement  combi 
nation  is  completed.  If  the  accuracy  measure  shows  otherwise,  take  another 
measurement  and  apply  the  variational  technique  to  calculate  the  wire 
currents.  Different  measurements  will  give  different  accuracy  measures.  Use 
the  measurement  which  gives  the  highest  accuracy  measure.  Continue  this 
process  until  a  satisfactory  accuracy  measure  is  reached.  Once  this  is 
reached,  it  is  essentially  the  best  measurement  combination  which  can  be  used 
for  future  stress  and  susceptibility  threshold  calculation  and  for 
determination  and  verification  of  LRU  specifications. 


VI.  CONCLUSIONS 


Analytical  procedures  have  been  developed  in  this  report  for  reducing  the 
number  of  measurements  required  for  a  reliable  aircraft  LRU  hardness  assess¬ 
ment.  A  deterministic  approach  via  a  variational  technique  has  been  under¬ 
taken  to  find  the  relationships  among  the  bulk  cable  current  and  the  individ¬ 
ual  wire  currents.  A  measure  of  accuracy  has  also  been  defined.  For  a  given 
accuracy  imposed  on  each  individual  wire  current,  optimal  measurement  ap¬ 
proaches  can  be  developed.  These  approaches  strongly  depend  on  the  excitation 
mechanism.  For  example,  under  two  simple  excitation  conditions,  the  optimal 
measurement  approaches  are  given  below: 

•  If  the  wires  in  a  cable  bundle  are  evenly  excited,  the  bulk  cable 
current  should  be  the  first  quantity  to  be  measured. 

•  If  the  excitation  of  one  individual  wire  is  dominant  over  the  others, 
the  current  on  that  specific  wire  should  be  the  first  quantity  to  be 
measured. 

From  the  measured  quantities,  the  individual  wire  currents  in  the  cable  bundle 
can  then  be  calculated  from  Equations  31  and  47. 

To  support  these  conclusions,  a  2-wire  cable  over  a  ground  plane  is  stud¬ 
ied  in  detail  to  find  the  effects  on  the  wire  currents  of  various  parameters 
such  as  the  excitation  mechanism,  bundle  configuration,  load  impedance,  source 
and  measurement  locations,  etc. 

The  results  derived  from  the  variational  technique  can  be  applied  to 
direct-drive  testing  to  determine  susceptibility  thresholds  and  to  develop  and 
verify  specifications  of  an  LRU.  For  example,  it  is  possible  that  a  proper 
excitation  of  one  wire  enables  the  determination  and  verification  of  the 
specifications  of  other  wires  in  the  same  bundle.  In  some  cases,  the  bulk 
cable  current  alone  can  be  considered  a  candidate  for  specification  of  an  LRU. 

New  experiments  have  been  suggested  for  verifying  and  improving  the  simple 
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APPENDIX  A 


RELATIONSHIPS  AMONG  WIRE  CURRENTS  WHEN  ONLY  BULK  CABLE  CURRENT 
AND  CURRENT  ON  ONE  STRONGLY  EXCITED  WIRE  ARE  MEASURED 


There  are  two  constraints  that  are  imposed  in  this  case.  One  is  the  sum 
of  individual  wire  currents,  Ig,  and  the  other  one  is  the  current  of  the  most 
strongly  excited  individual  wire,  1^.  Thus,  the  two  constraints  can  be 
written  as  follows: 


l,  >L  *  >B 

m=l  m 


\  =  h 

Lk  * 


Substituting 


I,  =  a  +  jb 

L  m  m 

m 

*B  *  *Br  +  jIBi 


Xk  =  l< r  +  jXki 


into  Equations  A1  and  A2,  obtains 


£  am  =  *Br 
m=l 

N 

I.  bm  '  hi 

m=l 


ak  =  !kr 


bk  '  "kl 


The  energy  function  of  the  Thevenin  equivalent  circuit  is  given  by 


N  N 


r  ■  I,  I,  iL  \  iL 

n=I  m*l  n  n.m  it 


r  ***-%! 


where  (Zt  )  =  (ZL  )  +  (Z$  )  .  Substituting  Equation  3A  into  Equation 

n,m  n,m  n,m 

10A,  gives 


p  ■  I,  l  ‘an  -  j6n'  Zt  (am  *  abm> 
n=l  m=l  n ,m 


(All) 


Now,  find  the  stationary  points  of  the  expression 

f  ■  I  I  <an  “  jbn}  Zt  {am  +  + 

n=l  m=l  n,m 


+  xl(IBr  "  am^  +  X2(IBi  '  J,  biJ 
m=l  m=l 


(A12) 


+  x3^kr  ‘  ak^  +  X4^ki  ‘  bk^ 


where  x1>  X2,  X3,  and  x4  are  Lagrange  multipliers.  Accordingly,  this  should 


result  in 


r—  =  0  ,  n  =  1,  2,  .  .  .  ,  N 
3an 

^  ’  °  ’  "  =  U  2 . " 


(A13) 


(A14) 


Consequently, 


I,  Zt  <a«,  *  Jb»>  *  l.  K  '  -  X1  '  x3 

11=1  n,m  m-1  m,n 


«nk  ■  0 


(A15) 


■S  l  h  <am  *  Jb»>  *  j  L  <an.  -  Jb»>  Zt  -  ‘2  *  Vnk  '  0 

m=l  n,m  m=l  m,n 


Assuming  that  (Zt  )  is  a  symmetric  matrix.  Equations  Alb  and  A16  as  can  be 
tn,m 

rewritten  as 
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*k,k 


Vm= 


=1  n=l 


m,n  J 


uk  ,m 


/ 


Substituting  Equations  A21  through  A24  into  Equations  A19  and  A20,  and  then 
using  Equation  A3,  one  finally  obtains 


APPENDIX  B 


ADDITIONAL  FIGURES 

This  appendix  presents  additional  figures  to  show  the  comparison  of  exact 
wire  currents  with  those  obtained  from  the  variational  technique,  and  also 
curves  to  demonstrate  the  effects  of  load  impedances,  source  types  and  loca¬ 
tions,  number  of  excited  wires,  line  configurations,  and  measurement  locations 
on  common-  and  differential -mode  currents. 

1.  COMPARISON  OF  EXACT  WIRE  CURRENTS  WITH  THOSE  OBTAINED  FROM  VARIATIONAL 

TECHNIOUES 

Figures  B1  through  Bll  illustrate  the  comparison  of  exact  wire  currents 
(Ij,  I2)  with  those  obtained  from  variational  techniques.  and  I£  are  the 
wire  currents  obtained  from  a  knowledge  of  IB  which  is  1^  +  Ig.  IJ  and  I£  are 
the  wire  currents  obtained  from  a  knowledge  of  one  individual  wire  current, 
which  is  taken  to  be  Ii . 


exp(-  at)sin(2n fQt)  mA 
a  «  5  x  106  <?_1 


.1 


I*(  tnfl  )  ,  Ij(  mfl ) 


VRRIRTI ONRL  PRINCIPLE  COMPARISON 


TIME  (  s  ) 

(b)  Time  history  of  1 2 >  1 2 •  and  *2* 

Figure  811.  Exact  wire  currents  (Ij,  I^)  and  estimated  wire  currents 
( 1 1 »  Ij;  Ij,  l'^)  using  two  variational  principle  techniques 
when  l  =  75  m,  x$  «  0.25  m,  xm  =  75  m,  Z^j  =  500  fl, 

=  250  Q.  and  =  100  Q. 


.  **  -r.  -j 


Z"*  ; 


2.  EFFECT  OF  LOAO  IMPEDANCE 

Figures  B12  through  B17  illustrate  effects  of  pin-to-pin  impedance 
(Z^^)  on  common-  and  differential -mode  currents  (Ig,  Iq). 
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EFFECT  OF  LORD  IMPEDANCE 


EFFECT  OF  LORD  I MPEDRNCE 
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(a)  Time  history  of  Ig. 
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30  m,  xs  =  0.25  m,  x^ 

Zm  =  2irfj(3xl0-7)  +  l/[2TTfj(7xlO'n)]f2  ,  and  lm 
2irfj(2.5xl0'7). 
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EFFECT  OF  LORD  IMPEDANCE 


Figure  B16.  common-  and  differential -mode  currents  ( Ig » for 

different  Z£12  when  l  =  30  m,  x$  =  0.25  m,  x^  =  30  m, 

Z. n  *  Z.„  =  100  +  2irfj(3xl0'7)  +  l/[2irfj  (7x1 0'1 1 )]  n 


TIME  (ps) 

(b)  Time  history  of  1^. 

Common-  and  differential -mode  currents  (Ig.Ip)  for 


different  Z£12  when  i  =  30  m,  x$  =  0.25  m,  x  =  30  m 
and  ZC11  =  Z12Z  =  2TTf j(3xT0‘7)  +  l/[2Trfj(7xl0"n)]  n 


IDC  mR)  .  IBCmR) 


<‘o  > 
n 


exp(-  at)sin(2itfQt)  V 
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a  =  5  x  10°  S-1 

f0  =  2  MHz 


mR  )  I»(  mR  ) 
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TIME  Cps) 

(b)  Time  history  of  Ip. 

Figure  821.  Common-  and  di fferential -mode  currents  for  different 

source  locations  when  £  =  in  m,  x  =  30  m,  Z„,,  =  Z  =  z 


IDC  mR)  ,  Ia(mR) 


In(  mR )  ,  IB(  mR ) 


(b)  Time  history  of  1^. 


gure  B25.  Common-  and  di fferential -mode  currents  (Ig.Ip)  f°r  different 
numbers  of  sources  when  i  =  30  m,  xs  =  0.25  m,  x,^  =  75  m, 

104  n,  Z,i9  *  «,  and  Z»,o  =  100  n. 


for  different 


IB<  mR  ) 


ID(  mfl  )  IB(  mR  ) 


Id(  mn  )  .  I»(  rnR ) 


2/2 


A0-A162  964 
UNCLASSIFIED 


A  VARIATIONAL  APPROACH  TO  DEVELOPING  RELATIONSHIPS 
AHONG  HIRE  CURRENTS  IN.  .  (U)  DIKEUOOD  SANTA  NONICA  CA 
N  ENGHETA  ET  AL.  NOV  G3  DC-FR-1G26.  S4G-1  RFML-TR-G4-97 
F29601-82-C-GG27  F/G  9/3 


END 


EFFECT  OF  LINE  CONF I GURBT I ON 


262.9  48.3 
48.3  317.9 


0.2  f 
(a)  Time  history  of  I 


262.9  48.3  \n 
48.3  317.9  /n 


1.4  0  .  S 

TIME  (  |j  s  ) 

(b)  Time  history  of  Ip. 


0 . 8 


Figure  B31.  Common-  and  differential -mode  currents  (Ig,I0)  for  different 

characteristic  impedances  when  t  *  30  m,  *  0.25  m,  xm  »  30  m, 

Z/H*  z/12  "  lpZ2  *  400  n»  and  U0  )  *  (Zc  ). 
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Figure  B32.  Common-  and  di fferential -mode  currents  (IB,IU)  for  different 
measurement  locations  when  l  *  30  m,xe  =  0.25  m,  ZM1  *  Z», 
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